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Abstract. - The adhesion of biological membranes is controlled by various types of receptor 
and ligand molecules. In this letter, we present a statistical-mechanical model for membranes 
that interact via receptor/ligand bonds of two different lengths. We show that the equilibrium 
phase behavior of the membranes is governed by an effective double-well potential. The depths 
of the two potential wells depend on the concentrations and binding energies of the receptors 
and ligands. The membranes are unbound for small, and bound for larger potential depths. 
In the bound state, the length mismatch of the receptor/ligand bonds can lead to lateral 
phase separation. We derive explicit scaling laws for the critical points of unbinding and phase 
separation, and determine the prefactors by comparison with Monte Carlo results. 



Introduction. - Biological membranes consist of a lipid bilayer with embedded or ad- 
sorbed proteins [1]. The adhesion of the membranes is typically mediated by various types 
of receptor and ligand proteins. Each type of receptor specifically binds to a complementary 
type of ligand in the apposing membrane. The adhesion of immune cell membranes is often 
mediated by receptor/ligand bonds of different length. The receptor/ligand complexes that 
mediate the adhesion of T cells, for example, have characteristic lengths of 15 or 40 nm [2]. 
During T cell adhesion, a lateral phase separation into domains that are either rich in short 
or long receptor/ligand bonds has been observed in the cell contact zone [3,4]. The domain 
formation is assumed to be driven by the length mismatch of the receptor/ligand bonds [5-11]. 

We consider here a statistical-mechanical model of two membranes interacting via long and 
short receptor/ligand bonds. In our model, the membranes are discretized into small patches 
that can contain single receptor or ligand molecules. The conformations of the membranes 
are described by the local separation of apposing membrane patches, and by the distribution 
of receptors and ligands in the membranes. We show that a summation over the receptor and 
ligand degrees of freedom in the partition function leads to an effective double- well potential. 
The potential well at small membrane separations reflects the interactions of the short re- 
ceptor/ligand bonds, and the potential well at larger separations the interactions of the long 
receptor/ligand bonds. The depths of the wells depend on the concentrations and binding 
energies of the receptors and ligands. 

We focus on the equilibrium phase behavior of the membranes, which exhibit two charac- 
teristic phase transitions. The first transition is the unbinding transition of the membranes, 
which is driven by an entropic membrane repulsion arising from thermal shape fluctuations. 
The second transition is lateral phase separation within the membranes, driven by the length 
mismatch of the receptor/ligand bonds. The length mismatch leads to a membrane-mediated 
repulsion between the two different receptor/ligand bonds, because the membranes have to 
be bent to compensate this mismatch, which costs elastic energy. This repulsion leads to a 
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Fig. 1 - (Left) A membrane containing long and short receptor molecules (top) adhering to a mem- 
brane with complementary ligands. In our model, the membranes are discretized into small patches 
(indicated by dashed lines). A membrane patch can contain a single receptor or ligand molecule. The 
membrane conformations are described by the local separation Zi of each pair i of apposing membrane 
patches. - (Right) Summing over all possible distributions of receptor and ligand molecules in the 
partition function leads to an effective double-well potential V E t for the membranes. The potential 
well at short separations z\ < Zi < Z2 reflects the interactions of the short receptor /ligand bonds, the 
well at larger separations zz < Zi < z^ reflects the interactions of the long receptor/ligand bonds. 



lateral phase separation for sufficiently large concentrations of the receptor/ligand bonds and, 
thus, sufficiently deep wells of the effective potential. 

For both types of phase transitions, we derive characteristic scaling laws. We confirm 
these scaling laws using Monte Carlo simulations, and obtain the numerical prefactors by a 
comparison with the simulation results. From the scaling arguments and simulations, we thus 
obtain relatively simple relations that characterize the phase behavior of the membranes. 

Model and effective potential. We consider two interacting membranes. One of the 
membranes contains long and short receptors, the other membrane contains complementary 
ligands (see Fig. 0). We use a lattice gas model and discretize the membranes into quadratic 
patches of size a x a [12]. Simulations with molecular membrane models indicate that the 
smallest bending deformation corresponds to a linear patch size a ~ 5 nm [13]. Each mem- 
brane patch can contain a single receptor or ligand molecule. We describe the distribution of 
receptors by a composition field ni with values 0, 1, or 2. The value ni = 1 indicates that a 
short receptor is present in patch i, ni = 2 indicates a long receptor, and ni = indicates that 
no receptor is present. The distribution of ligands in the apposing membrane is described by 
the composition field mi, with rrii = 1, 2, or indicating ligands of short receptors, ligands of 
long receptors, or no ligands in patch i of this membrane. The membrane conformations can 
be described by the local separation U of the apposing membrane patches i. 

The configurational energy of the membranes is the sum of elastic and interaction energies: 

H{z, n, to} = H c \{z} + H in t{z, n, to} , with H c i{z} = } \ (A d Zi) 2 (1) 

i 

To simplify the notation, we use the rescaled separation field Zi = (k/ 'a)\/ K/fcflT ', where 
k = k\K2/(ki + K2) is the effective bending rigidity of the membranes with rigidities n\ and 
K2 [14]. The discretized Laplacian A^Zi = za + Zi2 + z^ + zn — 4zi is proportional to the 
mean curvature of the separation field. Here, Zn to 2^4 are the membrane separations at the 
four nearest-neighbor patches of membrane patch i. The interaction energy has the form 

Hint = (Sn i ,lS mtyl Vi(z i ) + S nt: 2S„ H: 2V2(z i )~ 6 nit i(J,Ri - £ ni ,2/iR2 - 6 mj ,l Mil ~ Smt,2(J>L2) ( 2 ) 
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Fig. 2 - Monte Carlo data for membranes: (a) Critical potential depth U c of unbinding from the 
single- well potential (JSJ; (b) Critical depth for lateral phase separation in the symmetric double- 
well potential. The dashed lines represent the scaling laws HUH and 112H with numerical prefactors 
b — 0.025 ± 0.002 and c = 0.225 ± 0.02. The scaling laws are valid for large values of z\, z wc and Zba- 
The given value of b has been obtained from a fit to the Monte Carlo data points for z\ > 0.5 in (a), 
and the value of c is obtained from a fit to the data points for z wo > 0.5 and Zba > 1.4 in (b). 



where Vi(zi) is the interaction potential of a short receptor with its ligand, and V^Zj) is the 
interaction potential of long receptors. The Kronecker symbol o\j is equal to 1 for i — j, and 
equal to for i ^ j. The term 5 nit i5 mit i, for example, is only 1 for rii = 1 and m, = 1, i.e. when 
a short receptor and a complementary ligand are present in the two apposing membrane 
patches labeled by i. The chemical potentials of the short and long receptors are denoted by 
Uri and Hr2, and the chemical potentials of the ligands by fi^i and (J-L2- To simplify the 
notation, all energetic quantities, such as the configurational energy, interaction potentials 
and chemical potentials, are given in units of the thermal energy fc^T. 

The equilibrium behavior of the membranes is governed by the partition function Z. An 
exact summation over the degrees of freedom of the composition fields rij and m, in the 
partition function leads to 



Z = 



n 



IIE 



nE 



m =o J L i m i — o_ 



n 



-(Hel{*}+£ 4 V rf (*0) 



(3) 



with the effective potential 

V e f(zi) = - In [Co + e"" 1+ ^ L1 (e- Vl ^ - 1) + e ^+Mi.2 ( e -%(«) _ x )l (4) 



and Co = (1 



0(1 



^ L2 ). Here, we characterize the molecular interactions 



by the square-well potentials V\(zi) = —U\ for z\ < zt < Z2 and otherwise, and Vi^zC) = — C/ 2 
for Z3 < Zi < Z4 and otherwise. For zi < Z3, the effective potential V c f(zi) is the double-well 
potential shown in Fig. ^ with the effective depths 



Uf = In 



l)/Co 



and If? = In 



1 



(e^-l)/Co (5) 



of the two potentials wells with respect to V e f(zi = 00) = — In Co- 
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Unbinding from single-well potential. - We first consider the single-well potential 

V(z) = — U for zi < z < Z2 and V(z) = otherwise (6) 

A potential of this form is obtained by setting Uf or J7f equal to for the double-well 
potential shown in Fig. ^ The critical scaling behavior of membranes interacting via such 
a potential is similar to the scaling behavior of strings. Since strings can be studied with 
analytical methods, we first determine the critical potential depth U c of strings as a function 
of zi and Z2 and subsequently compare with Monte Carlo results for membranes. 

Strings are lines governed by tension. We consider here two interacting strings in two 
dimensions. The conformations of the strings can be described by the local separation l{x) 
perpendicular to a reference line. The strings are, on average, parallel to this line. To simplify 
the notation, we use again a rescaled separation field z(x) = l(x)y/ 'a /ksT where a is the 
effective tension of the two strings. The effective Hamiltonian has the form [15] 



H {z} = J dx 



(7) 



where V(z) is the interaction potential of the strings. In the continuum limit, the separation of 
two interacting strings is equivalent to the spatial coordinate of a quantum-mechanical particle 
in one dimension [15]. The free energy of the strings then corresponds to the ground-state 
eigenvalue Eq of the Schrodinger-type equation 

-^T+V{z)fl> k {z)=E k fl> k (z) (8) 
oz z 

In the case of the square- well potential © , the eigenfunction ipo has the form 

ipo(z) — Ai(exp(fcz) — exp(— kz)) for < z < z\ 

= Ai cos(az) + A3 sin(az) for z\ < z < z^ (9) 
= A4 exp(— kz) for z > zi 

with a 2 = Eq + U and k 2 = —Eq. Since the strings cannot penetrate each other, we have 
tpo(z) = for z < 0. The two strings are bound for E < 0, and unbind at E — 0. The 
eigenfunction ipo(z) and its derivative <hpo(z)/dz have to be continuous at z = Z\ and z = Z2. 
These four continuity conditions lead to a transcendental equation. At the unbinding point, 
i.e. at Eq = 0, the transcendental equation is cos [(Z2 — zi)y/Uc\ = ziy/TJcSin [(z 2 — Z\)^/Uc\ . 
For zi = 0, this equation simplifies to cos(z2\fTT c ) — 0, which implies U c — 7r 2 /(4z|). For 
zi 7^ 0, the equation can be rearranged to Ziy/TJ^t&n Uz% — zi)y/U c \ — 1, which leads to 

Uc = —r-^ r for (Z2 -zi)y/U c *L\ (10) 

Z 1 (Z2 - Zl) 

with prefactor 6=1. 

Functional renormalization arguments indicate that membranes have the same critical 
scaling behavior as strings [14]. It is therefore reasonable to assume that the critical potential 
depth U c of two membranes interacting via the single- well potential © is governed by the same 
scaling law H10|l . but with a prefactor b that is different from the factor 6=1 for strings. To test 
this proposition, we determine the critical potential depth for membranes with Monte Carlo 
simulations. In the simulations, we use the discretized Hamiltonian 7i{z} = H c \{z} +Y1 V(zi) 
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with the elastic energy (JIJ and the single- well potential JJjJ, and attempt local moves in which 
the separation Zi of patch i is shifted to a new value Zi + C where Q is a random number between 
— 1 and 1. Following the standard Metropolis criterion [17], a local move is always accepted 
if the change ATL in conformational energy is negative, and accepted with the probability 
exp(— A7i) for AH > 0. We perform simulations with up to 10 7 attempted local moves 
per site i and membrane sizes of TV = 120 x 120 patches. At the unbinding point, the 
correlation length and the autocorrelation time of the membranes diverge towards infinity. To 
determine the critical potential depth U c , we perform simulations for larger potential depths 
U > U c at which the lateral correlation length is significantly smaller than the membrane size. 
Thermodynamic averages (f ) and (Pf>) for the mean separation z = 5Zi=i z i anc ^ the fraction 
p, of membrane patches bound in the potential well then do not depend on the finite system 
size. The critical potential depth U c is obtained from extrapolating (z) and (P,) as functions 
of U to the critical values l/(z) = and (P&) = 0. Monte Carlo results for U c at various values 
of the potential lengths z\ and zi are shown in Fig.[5Ja). Since the membranes are discretized, 
we only expect agreement with eq. (11011 in the continuum limit of large z\ and Z2- From fitting 
eq. ifirljl to the Monte Carlo results for large z\ and 22, we obtain b = 0.025 ± 0.002. 

Lateral phase separation in a symmetric double-well potential. - We now consider the case 
where the two wells of the effective potential shown in Fig. ^ are located at large membrane 
separations and neglect the hard- wall repulsion at z = 0. By comparing with the results of 
the previous section, we will later show that this case is realistic for the typical dimensions of 
cell receptor/ligand bonds. We assume here that the two potential wells have the same width 
z wc . For symmetry reasons, the lateral phase separation then occurs at the same depth of 
the two potential wells. Therefore, we focus here on the symmetric double-well potential with 
U cf = Uf = U^- Our goal is to determine the critical potential depth Ujf for lateral phase 
separation as a function of the width z wc and separation Zba of the potential wells. 

To derive a scaling relation for Ufr , we first show that the free energy of membranes bound 
in a single-well potential of depth U and width z we scales as 

T ~ {Uz wc f (11) 

for small values of U and negligible hard-wall repulsion. This scaling law can be obtained again 
from the analogy with strings. If the hard wall is negliglible, the ground-state eigenfunction i/jq 
is symmetric with respect to the center z m — (z2 + z\)/2 of the single well. The eigenfunction 
given in eq. © then simplifies to ipo{z) = B\ exp(— k\z — z m \) for z < z\ or z > Z2, and ipo(z) = 
B-2 cos(a(z — z m )) for z\ < z < Z2- From the continuity conditions at z\ and Z2, we obtain the 
transcendental equation ^E + U tan (i\/Po + Uz we ) = y/—E . For small values U, the free 
energy Eq of the strings is also small. We thus can expand the tangent in the transcendental 
equation and obtain the explicit solution E = (— 2 — Uz^ c + 2y/l + Uz^ e )/z^ e ~ — (t/z we ) 2 /4 
to leading order in U, and, thus, the scaling law for the free energy T = Eq as in (|ll|l . 

We now come back to the symmetric double-well potential with well depth U ci . According 
to ljll(l . the free energy of the state where the membranes are bound in one of the potential 
wells scales as T ~ (U cf z wc ) 2 . More precisely, T is the free energy difference with respect 
to the unbound state. This free energy difference includes the entropy loss of the bound 
membranes, and is therefore different from the binding energy U ci for the wells. The absolute 
value, \J-\, corresponds to a free energy barrier that the membranes have to cross from one 
well to the other well. A central, previous result is that the membranes exhibit a lateral phase 
separation if this free energy barrier exceeds a critical barrier height T c ~ l/zj~ a [16]. From 
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Fig. 3 - Phase diagram of membranes adhering via long and short receptor/ligand bonds. The 
membranes are unbound for small well depths and ?7| f of the effective interaction potential shown 
in Fig. Q i.e. for small concentrations or binding energies of receptors and ligands (see eq. 0). At 
large values of and £/| f , the membranes are either bound in well 1 or well 2, i.e. they are either 
bound by the short or by the long receptor/ligand bonds. At intermediate well depths Uf and [7f f , 
the membranes are bound in both potential wells. The critical point for the lateral phase separation 
(star) follows from eq. (II 12 P - For typical dimensions of cell receptors and ligands, the critical well 
depth Uc C for lateral phase separation is significantly larger than the critical depths of unbinding (see 
text). The critical unbinding points for Uf 1 — or J7f f = (dots) follow from eq. I1U1 , 



this scaling relation for the critical barrier height and eq. we obtain the scaling law 

Uf=^— (12) 

ZweZha 

for the critical depth of the double-well potential. 

To verify eq. 1)12(1 and to determine the numerical prefactor c, we compare again with 
Monte Carlo results. The critical potential depth can be determined from the moments C2 = 
(z 2 )/(\z\) 2 , and C4 = (z 4 )/(z 2 ) 2 . Here, z — jjJ2iLi z i ^ s ^ ne spatially averaged separation, 
and (• • • ) denotes an average over all membrane configurations. The values of these moments 
depend on the correlation length £ and the linear size L of a given membrane segment. At 
the critical point however, the correlation length £ diverges, and the values of the moments 
become independent of L [17]. The critical depth Uf of the symmetric double- well potential 
can be estimated from the common intersection points of the functions C2(U cl ) and C4(U ei ) at 
different values of L [17-19]. We have considered quadratic membrane segments with a linear 
size of L = 10, 20, and 40 patches and periodic boundary conditions. From the comparison 
with the simulation results shown in Fig. H£b), we obtain c = 0.225 ± 0.02. 

In the present study, we ignored possible effects of membrane tension, a, which dominates 
over the bending energy on length scales larger than the crossover length \J k/<j, but is neg- 
ligible on smaller length scales. For lipid bilayers with a bending rigidity k of about 20 fcsT 
and tensions of a few /zJ/m [20], the crossover length \fnja attains values of several hundred 
nanometers, which is significantly larger than the membrane deformations between domains of 
short and long cell receptor/ligand bonds [11]. Therefore, tensions up to a few //J/m should 
only weakly affect the critical potential depth l|12|) for phase separation. However, tension 
prevents unbinding from the short-range potentials considered here [12]. In addition, active 
switching processes of stickers can have a strong impact on the unbinding behavior [21]. 
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Phase diagram and Conclusions. In the case of T cells, the long and short recep- 
tor/ligand bonds have linear extensions of I = 40 and 15 nm, respectively [2]. Reasonable 
values for the width of the potential wells range from l wc — 1 to 5 nm. The rescaled width z wc 
of the two potential wells is then significantly smaller than rescaled separation z\ of the first 
well from the hard wall and the separation Zba of the two wells, and the ratio of z\ and Zba 
is around 0.6. For these values, the critical potential depth llf for lateral phase separation is 
much larger than the critical depths of unbinding. According to eq. (JTDJ, the critical depth 
of unbinding from well 1 is U c i — bfz\z wc for U% f = 0, and the critical depth of well 2 is 
U c ,2 = b/z 3 z we for Uf = 0. The critical depth U° f = c/z^ a z wc for lateral phase separation 
is much larger than both U c \ and f/ c 2 since the ratio of the prefactors in the critical scaling 
laws <|10|1 and (|12|l is b/c ~ 0.11. For Uf — C/| f = U]~, the membranes therefore are clearly 
bound in both wells, which justifies that we have neglected the hard- wall repulsion in deriving 
eq. <|12|) . The resulting phase diagram is shown in Fig. EI 

In summary, we have characterized the phase behavior of membranes adhering via long and 
short receptor/ligand bonds. We have mapped the partition function of such membranes on 
the corresponding partition function of homogeneous membranes interacting via an effective 
double-well potential, and have determined the critical points of unbinding and lateral phase 
separation as a function of the characteristic potential parameters. These parameters are the 
effective depths of the two wells, which depend on the binding energies and concentrations of 
the receptors and ligands, and the characteristic length scales of the potential, which depend 
on the dimensions of the receptor/ligand bonds. 
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